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Abstract: 

The classical and quantum aspects of planar Coulomb interactions have been studied 
in detail. In the classical scenario, Action Angle Variables are introduced to handle rela- 
tivistic corrections, in the scheme of time-independent perturbation theory. Complications 
arising due to the logarithmic nature of the potential are pointed out. 

In the quantum case, harmonic oscillator approximations are considered and effects of 
the perturbations on the excited (oscillator) states have been analysed. 

In both the above cases, the known 3+1-dimensional analysis is carried through side 
by side, for a comparison with the 2+1-dimensional (planar) results. 
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I. Introduction 



In recent years, physics of the lower dimensional systems are playing more and more 
important roles, not only as toy models for the physical world, but also as an indepen- 
dent entity. On the one hand, effectively 1+1 and 2+1-dimensional systems have been 
constructed, while on the other hand, rich theoretical structures have been unearthed in 
these low dimensions. 

In the present work, we have studied in some details, the classical and quantum aspects 
of planar Coulomb interactions. We have also attempted to estimate effects of relativistic 
corrections. The Darwin Lagrangian (and the Hamiltonian) , that we have considered 
has appeared before in 0. Their roles have been further elaborated in ||, where the 
Coulomb interaction between anyons has been considered. These studies are done in 
the spinning particle model of anyon, proposed by us ||. Unlike the 3+1-dimensional 
Coulomb problem ||, the logarithmic nature of the Coulomb potential poses computa- 
tional as well as conceptual problems. Instead of trying for a full solution of the problem, 
we obtain the (degenerate) frequency in the Action Angle Variable (AAV) formulation |6j , 
for a restricted set of orbits. Then we apply the time independent perturbation theory 
in the lowest order in rj — v/c, (v and c being the velocities of the particle and light), 
for heavy and slowly moving particles, to study the relativistic corrections. Once again 
these steps are not straightforward and we have to make some ansatz which are justified 
numerically in a later chapter. We have shown that the perturbation lifts the degeneracy. 
However, the relations connecting the original coordinates to the AAVs are not well de- 
fined for the whole range of variables. Hopefully, this problem might turn out to be more 
of a computational nature. 

In the quantum case, we have started by performing a harmonic oscillator approxi- 
mation about the potential minimum and have compared the 3+1 and 2+1-dimensional 
results. Interestingly, the binding energy obtained from classical and semiclassical (or 
quantum) approaches are not widely different in the planar case, whereas they differ by 
orders of magnitude in the 3+1-dimensional case. The classical and quantum cases differ 
in the value of the angular momentum of the particle. In the former, we obtain it from the 
particular potential profile we have considered, whereas in the latter, we simply assume 
it to be of the order of h. Next we study the effects of perturbation in the harmonic 
oscillator excited states JF], which contribute only from the second order in rj. Here we 
face a conceptual problem concerning the range of the radial coordinate variable. The 
correction depends on the boundary condition, where we fix the zero of the logarithmic 
potential. 

As the foreward indicates, it may seem that our analysis has both solved and unearthed 
problems, (at least as far as the study of perturbations are concerned), of equal number! 
This obviously underlines the fact that we must know the classical and quantum aspects 
of planar potential problems more thoroughly. Also it is worthwhile to mention that some 
of the troublesome questions are inherently present in 2+1-dimensional system, which are 
not a reduction of an original 3+1-dimensional system, (with a translation symmetry in 
one direction), as is normally encountered. 

The paper is organised as follows. In section II, we reproduce the 3+1-dimensional 
(i) Coulomb problem results in terms of AAV. This actually outlines the approach to be 
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pursued for the planar Coulomb potential, which is carried out in section III. The 

validity of the ansatz, introduced in section III, are checked in section IV. Section V 
constitutes a brief analysis of the time independent perturbation analysis of the relativis- 
tic corrections of the present model, where some problems in the numerical study have 
been stressed. Section V comprises of the quantum analysis where again the logarithmic 
potential raises some non-trivial questions. The 3+1-dimensional analysis is also carried 
out side by side for completeness. 



II. Classical Analysis in 3+1-dimensions 



In this section, we briefly discuss the text book problem || of 3+1-dimensional Coulomb 
problem in (AAV), to outline the approach and for comparison with the 2+1-dimensional 
case at various stages. 

For the cases where explicit solutions of the time evolution of coordinates or the equa- 
tions of orbits are not required, (or are not computable), time independent perturbation 
theory comes to the rescue. If one can express the original coordinates and momenta in 
terms of AAV, the perturbations in the frequencies (of periodic systems) are computable 
in a systematic way and with little fuss. However, the catch is in the former task! As 
is well known this poses no problem in 3+1- dimensions, but the analogus problem runs 
in to severe computational trouble in 2+1-dimensions owing to the logarithmic Coulomb 
potential. 

The Lagrangian and Hamiltonian of two interacting particles of mass 2m each and 
charge Q and — Q are, 

T 1 I -2 i 2 /' 2\ i Q {-i \ 

- -m(r + r ip ) + 



2 v Y ' 47re r' v ' 

p 2 I 2 Q 2 

2m ^2mr 2 47reor + e //' 

where the r and ip are the relative coordinates in the plane of the motion and m is the 
reduced mass. The conjugate momenta are 

p r = mf, = mr 2 ip = I = constant. 

In Figure (1), ^rV e ff — — ^) vs. r is plotted, where A = metre = 0.01 metre, 

for later comparison with 2+1-dimensional results. 

The action variables in a general r, 6, <p coordinate system are 

■h = <j> Pkdq k} k = r, 9, 0, (3) 

the integral being a line integral over complete periods of orbits in (q^, p^) plane. The 
generic central force Hamiltonian is 

H = ^{pl + - 2 +~^) + V{rl (4) 
2m r z r z stn z t) 

where pk = Jjp. The total angular momentum I of (|2|) is related to this notation by 

p 2 

I 2 = pi -\ — r ~^7a • Here one obtains, 

Jj, = 27rp^, J e = 2ix{l- ptj,), (5) 



3 



J r = i Pr dr = 2 [ r \2mE + ^ - { -^f)^dr, (6) 

where, r\ and r2 being the two roots of the integrand, are the turning points of the classical 
motion and V(r) is the Coulomb potential. E is the total energy of the system. This 
integral can be done directly or via contour method, with the result, 



= ~(J, + J*) + j^zg. (7) 

Now (0) can be inverted to write in terms of the action variables only, 

mQ A 



E 



8e 2 (J r + J e + Jj>y 



The degenerate frequency v k = jjj- is, 



By performing the indefinite integrals in 



W 



r(j) rO rT 

I p^dcj) + / p e d6 + / p r dr, (10) 



one can see that Hamilton's characteristic function is completely separable. The angle 
variables are defined by Wk = §j- and inverting these relations the original coordinates 
are expressible in terms of the AAVs. Note that all the integrals are trivial here. From 
here the road to time independent perturbation theory is clear ||. The procedure will be 
elaborated in the main body of our work. 

On the other hand, one can get the crucial result of in an easier way, albeit in a 
more restricted setting, by the following argument. The energy of the circular orbit is, 

E C = H \ Pr=0 , _ = - 8eg(Je + ^ )2 , (11) 

where r c is the radius of the orbit, r c = 2A, obtained from the condition — | r= r = 0. 
For circular orbit we have p r = and hence J r = 0. 

Now if one is assured of the existence of closed orbits, (for example by invoking 
Bertrand's theorem 0), the frequencies must be the same or at most rational multiples 
of each other. Hence, at least for the completely degenerate case, for non-circular orbits, J r 
must enter the energy expression in the same way as Jq and in (|TTD and so one recovers 
the general result (|8|), by simply replacing ( J e + J^) by J r + Jg + in ([TT|). Obviously the 
former rigorous result is much stronger, showing that there are no non-degenerate orbits 
in the 3+1 dimensional Coulomb problem. We will try to apply the latter scheme in the 
2+1-dimensional problem under study. 

III. Classical Analysis in 2+1-dimensions 
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Let us start with the Coulomb Lagrangian and Hamiltonian, 

C=^(P + rW)-^ln^, (12) 
2 2ne r 

n = f~ + + In-) = KE+ V eff , (13) 

2m 2ireo Q mr z ro 

with pq = mr 2 8 — I, a constant, and p r = mr. The Coulomb potential is consistent 
with the Gauss law, V.E = ^, where p is the charge density, E the electric field and eo 
(or "permittivity") is just a property of the surrounding vacuum, which is related to the 
other property /iq, ("permeability"), by the relation eo/io = Note that now e is such 
that ^- has dimension of energy. The Coulomb potential is attractive for r < ro, ro being 
the distance scale where the potential vanishes. 

In Figure (2), we have shown the potential profile of V e ff = f° r B = T^pr- 

The graph is obtained by plotting ^w-V e ff vs. r, with r = 1 and the value of the 
parameter B = 1CT 5 (metre) 2 . Note that B is structurally similar but dimensionally 
different from the analogous parameter A in 3+1-dimensions. 

The action variables are, 

Je = j> Ped6 = 2tcI, J r = <j> p r dr . (14) 

The straightforward computation of J r as in @ is not possible since the transcendental 
nature of the integrand forbids explicit evaluation of the roots (or turning points). A 
contour integration might help but we have not been able to perform it. Instead we take 
recourse to the qualitative analysis described at the end of section (II). 
The energy and radius of the circular orbit are, 

eo Je 



2B . (16) 



2ixm Q 

Bertrand's theorem || assure us that closed orbits are allowed. Then at least for the 
degenerate orbits, where u r — i/q, the energy is, 

E = -^-(i + ln[ s H°~ Je + Jr ]) . (17) 
27re V 2 [ V 2nm Qr U V ; 

By inverting this relation we get, 

r ,2%e E 
—exp{ — 
a Q 



J 9 + J r = —exp( ^ 2 ) 



where a = y 2nmQ 2 an< ^ 6 * s ^ ne exponential number. We will give numerical estimates of 
the action integral in ( [L8|) later to show the validity of the above ansatz. Hence the single 
non-trivial frequency is, 

Vr = Ve= 7T—{ j i t ) = 7T—T , (19) 

27ie J r + Je 271€qJ 2 
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where we have defined a new set of action variables, Jq = J\ and J r + Jg = J 2 . The above 
is the unperturbed degenerate frequency. 

The next task is to obtain expressions for the coordinates r and 6 in terms of AAVs. 
For this we define the Hamilton's characteristic function, 



/9 rr 
p e d6 + / p r dr 



4+^-74-^. (20) 
2n J 4-K z r z neo r 

Note that the arbitrary ro cancells out in the right hand side, but we carry it to keep an 
account of the dimensions and finally take tq = lmetre. 
The two angle variables are, 



i 



dW 6 Ji r dr , ^ J? Q 2 m , r . 
oJi 2tt 2tt z J r z 4n z r z neo r o 

dW Q 2 m r , , n ^ J? Q 2 m , r , i 
w 2 = — = ^— / dr{2mE - - 2 — In- )"s 22 

Once again we are stuck with integrals having logarithm inside. Let us rewrite the di- 
mensionless angle variables in the form, 



9 HT r dr r B 



w i = — ~ °Ji\l— I "a (« - In- - -j) ^ , (23) 



1 \[2~e f , , , r B 

w 2 



2n aJ 2 J r 
where 

aJ 2 eoJx 2 (aJi) 2 

a = In , B = — = 

ro AnmQ 2 2e 

We use the following arguments to make the substitution, 



/ dr(a - In -)s , (24) 

J rn r 2 



a - In— - ^ = (r - n)(r 2 - r)^-^n 2 , (25) 
r cr J 2 



where r\ and r 2 are the smaller and larger roots and n is an integer, to be determined 
later. From the nature of the potential well in Figure (2), we know that for the energy 
range > E > E c , there are only two roots, (or turning points), of the expression of p r in 
(P0|). Since in the classical regime, we have to stay inside the turning points for physical 
motion, the replacement in (|25|) is justified. Later we will compare the two expressions 
in ( |25|) numerically. The constant factors are picked in such a way that when (|2"0f) is 
inverted, the libration coordinate r becomes a periodic function of w 2 . From (p4[) and 
(p5[ ), we find, 

1 /" e?r , x 

W2 = . (26) 

zrm J /( r _ n)(r 2 — r) 

Next we come to the roots ri and r 2 . Indeed, numerical values of these roots will not 
suffice, since we require analytic expressions in terms of AAVs. 



6 



Computing the larger root r 2 is easier, in the limit r — > oo, we get 

r 2 = r exp(a) = oJ 2 . (27) 

However, estimation of the smaller root 7*1 is tricky. We use the algebraic consistency 
of the above relations, (in the special case of the parameter value we have chosen from 
Figure (2)). From fll8|) and ([25D, we have, 



y/e 2nmQ 2 p r B * 

J r = 1 \\ / ar(a — In 

V27T V ee Jri r H 

2ne f 1 " 2 r — r en 2 

— — / dry/(r -n)(r 2 - r) = — — -(r 2 - n) . (28) 



Using ( p7|) we simplify ( p8|) to get, 

Jr =^ W + i-^) = ^ (J - + J 1 + ^-^). (29) 
4 a z J 2 cr 4 a z J 2 a 

Hence comparing left and right hand sides, we must have, 

— = 1, r 2 1 -2ar 1 J 2 + (T 2 J 1 J 2 = . 
Solving for the nearest integer value of n and for r\ we get, 

n^2, n = aJ 2 ±aJ 2 (l-^ . (30) 

When we give the numerical estimates in the next section, it will be seen that J 2 > J\. 
Since aJ 2 — r 2 > r±, we find 

n~4 + sf). (3D 

Later we will exihibit values of the roots obtained numerically and obtained from our 
approximate analytical expressions, for a comparison. 

Taking all this into account, we finally get the angle variable as, (with n = 2), 

1 . _j2r- (ri + r 2 ) , , 

w 2 = —sin , (32) 

47T T 2 — T\ 

9 ct 2 JiJ 2 f dr 
wi = — 



(33) 



2n 2\/2e J r 2 ^J(r — r 1 )(r 2 — r) 

9 a 2 JiJ 2 [-rir 2 + (ri + r 2 )r-r 2 ]2 
2n 2\/2e rr l r 2 
By inverting these relations we derive the desired expressions, with 47ru> 2 = 4>, 

r = ^[(r 2 - ri)szn<j) + (r 1 + r 2 )] 

= y[(l-^W+(l + P > (34) 
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6 1 / 1 — sirup , 

— = w\ H ^-W t rr^(10 + Ssinm + 2sm 6 — 3a + asino) , (35) 

where p = is a small parameter. Before applying these relations in perturbation theory, 
, in the next section we will demonstrate the validity of our assumptions numerically. 

IV. Validity of the Ansatz - A Numerical Study 



In the previous section, two crucial assumptions were introduced, in a rather cavalier 
fashion: 

(i) The energy expression of for a circular orbit in fllTf ) was generalized to any degenerate 
closed orbit in (|7]). 

(ii) The expression a — — ^ was replaced by a quadratic polynomial in r in (|2T 
Below we give the numerical comparison by going through the following steps. 
(A): From Figure (2), we obtain a reasonable value of the parameter 



For this particular B, 



B 



E P . 



AnmQ 2 



10" 



'IB 



Q 1 



+ ln(VW)) 



Q 2 



-y 



27reo 2 27reo 
We have put r = 1 metre. Hence the range of E, such that only two turning points 
occur, is > E > E c = JSL-(-4.96). In Table (1) below, we choose some typical values 
for E, (where 2n J$ E = y), and the corresponding roots r\ and r 2 from Figure (2). 



y 


n 


r 2 


s/Yej 1 


-i 


0.00133 


0.3678 


0.261 


-1.5 


0.0014 


0.22308 


0.154 


-2 


0.00149 


0.13523 


0.089 


-2.5 


0.00159 


0.08196 


0.051 


-3 


0.00172 


0.04958 


0.027 


-3.5 


0.0019 


0.02985 


0.013 


-4 


0.00216 


0.01773 


0.006 


-4.5 


0.00264 


0.01006 


0.001 



(B): First of all, let us check the smallness of p = j- 

E = and E c = £^ln(V2eB), p = V2^B w 0.0074 and p 
J\ — Ji. We have, 

x 



UnBmQ 2 V2^B 



For the two extreme values of 
1, (since for circular orbits, 

0.0074 



Jo 



ee 



exp( 



a 

2neoE, 
Q 2 ' 



a 



a 



1 ( 2ne E 
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We now compare the above results with our analytical expressions for r\ and r 2 in (|3T| 
and in fl27l) . Let us take a generic value E a and the roots rf 2 . Thus we have E a = 
We also have 

J i = -> ^2 = —exp( — — — ) = -exp(y) p = xe y . 
a a Q z a 

Hence the roots are, 

aJi aJ\ x x 2 e~ y 

ri = 1 a = — I , 

2 8 2 8 

r 2 = aJ 2 = exp(y) . 

This are listed in Table (2) below. 



y 


n = f(i + =p) 


r 2 = e y 


2CJT2 


-l 


0.00371 


0.36788 


0.49 


-1.5 


0.00372 


0.22313 


0.293 


-2 


0.00374 


0.13534 


0.174 


-2.5 


0.00377 


0.08208 


0.102 


-3 


0.00382 


0.04979 


0.058 


-3.5 


0.00391 


0.0302 


0.033 


-4 


0.00406 


0.01832 


0.015 


-4.5 


0.0043 


0.01111 


0.006 



(C): We here show the comparison between the two forms of J r . One is obtained by the 
following numerical integration, with the limits (obtained numerically) in Table (1). 



lQ 2 m 
ne 



r2 dr( 27ie ° E - neoJ ^ 
i Q 2 4:Ti 2 mQ 2 r 2 



In—) 2 



lQ 2 m 



TCt Jri 



Q 2 



B 



In- 



r 



2e 



On the other hand, we have the analytical result from (ESI), 



en 
Aa 2 J 7 



(r 2 - nj 



e(r 2 - n) 
2r 2 a 



(36) 



(37) 



This are displayed in Table 2, in order to compare between \ n and e ^ r2 2 ~ 2 ri ^ for 

different values of E and corresponding r% and r 2 , from the two Tables. 
(D) In Figure (3), comparative studies of the substitution of the logarithmic function of 
r by the quadratic polynomial in r , in (|25|), is shown for different values of the energy. 
In fig. (3) we plot the two functions in ( p5|) inside the range r" to r% for four values of E a , 



f(r) 



a 



, r B 
In 

ro r A 



In 



a Jo 



,aJ\. 
1 2er ' 



y — ln(r) 



x 



4e 2 r 2 
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Note that even though the natures of the curves are somewhat different, (the respective 
peaks being in different places), the results are reasonably satisfactory, since we only use 
the function inside the action integral and effects coming from the spatially separated 
maxima average out. 

V. Time Independent Perturbation Theory 



In this section, we only reproduce the expressions required for first order perturbation 
of the frequency ||, where i] = - is our small parameter. The perturbation term Hi in 
H = H + rjHi is expressed in terms of the unperturbed AAVs, 

ai(J ) = Hi(w , J ) , (38) 

where overbar denotes averaging over a complete period of Wo, with v = w = 
uot + Bo, Bo being a time independent constant. After the averaging, it is legitimate to 
replace Jo by J, (which is correct to 0(rj)), and one derives the frequency u, correct upto 

o(v), 

da da x 
U= dJ = U0 + r] ^J- (39) 
The perturbation terms that we consider has been discussed in detail in and [[| . 

H = Ho + V Hi = [- + lr-ln-\ + [a(l + In-) + -^-} , (40) 
m ivreo r r mcr 

where a = and b = 6 Jr%Pj are dimensionless variables and j is the spin of the free 

mcr mcr J 1 

particle. In it has been shown that in interacting systems of this form, there is a 
screening effect on the particle spin. The higher order terms 0] in a and b are dropped. 
We replace Pi by mvi and note that both a and b are of 0(rj). It should be kept in mind 
that the a dependent term (and higher order ones present in the original work 0) are 
conventional ones but the j dependent b terms are the spin effects. Here our aim is to 
check the stability of the closed orbits once the perturbation is switched on. 

It is immedietly apparant from (|34|), (|35| ) and fl4"0| ) that the degeneracy in the frequency 
is lifted. But we also notice that due to the necessary averaging to be done with respect 
to the angle variable (ft = 47ru> 2 , the expression for 6 in ( |5oD is not well defined, as the 
denominator may vanish for some value of the angle. It should be pointed out that the 
above mentioned problem crops up even in the j independent correction terms as well. 
Hence it is not clear whether the fault lies in the general perturbation scheme or in the 
several approximations that we have introduced. 

VI. Quantum Mechanics 



We start this section by giving estimates of the binding energy of a quantum system in 
3+1-dimensions. Following the analysis in ||, the total energy expression to be minimised 
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IS 2 9 

E( r = R) = - , (41) 

v ; 2mR 2 4ne R V ' 

ehere R is the ground state "radius" of the atom, h is a typical value of the angular 
momentum in the quantum domain. This gives us 

47re ft 2 _ 87r 2 h 2 _ 2 

il = -rr— = <T , 



E{R) 



mQ 2 

mQ 4 _ Q 2 e 



2.(47re ft) 2 647r 3 e ft 2 a 2 ' 

This comes to roughly R « 0.5291CT 10 metre and E w —13.5 eV. The exactness of the 
numerical value of E is merely a coincidence ||. 

Following this approach in 2+1-dimensions, we minimise, 

E( r = R) = -0— + -^-In-. 
y ' 2mR 2 27re r 

From this we derive, 

R = —=a, E = -*—ln . (42) 

Ve 27re r 

The parameter a is expressionwise identical in both cases but differ in dimension and is 
the introduced in section (III). 

It is amusing to note the following point: From Figure (1) and Figure (2), we can get 
an order of magnitude idea of the angular momenta involved; 

y/eA 

I = w 2.1CT 30 Joule — Second, 3 + 1 — dim., 

2ixo 

e\f2B 

I = — . 2 + 1 - dim. (43) 

2ixo 

In estimating the numerical value of I in 3+1-dimension, the particle is taken to be an 
electron. Note that (at least in 3+1-dimension) these values of I are much larger than 
h, the typical order of quantum angular momentum. Obviously this I value in 3+1 
dimensions will drastically change the ground state energy from the correct one, making 
it much smaller. However, since the angular momentum occures in the logarithm in the 
2+1 dimensional energy expression, there is no order of magnitude change in the energy 
value. 

One can do a harmonic oscillator approximation about the potential minima in Figures 
(1) and (2). We will follow this up in the 2+1 dimensional case with effects of the 
perturbations introduced in the previous section. 



In 3+1 dimension, we expand the effective potential energy about r c 



toff 



)2 ! 



Q 2 ( 1 | e J 2 1_ 

47re r c + x 2nmQ 2 (r c + x) 



Veff = ^(-r4n + hvr ) . (44) 
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and identify the quadratic x-term with the harmonic oscillator potential, and obtain the 
frequency , 

Q 2 eirQ 2 1 irmQ 4 

\Ane mrl 47ie Jf a 2 2e 2 Jf ' 1 ' 

Hence the excitation levels above the ground state are 

E n = nhQ , n — 1, 2... . 
Identical analysis in 2+1 dimension leads to an expansion of V e ff about r c — 

Q 2 r c + x e Jl 

Ufl -|- 

27reo r 27imQ 2 (r c + x) 

which produces the frequecny 



nm Q ' 



V eff = + T^ZT^TZ I ' ( 46 ) 



Q = t = ^ 2 (47) 

2 v / 27rmJ 1 cr 2 v^eoJi 

Note that Q is independent of the particle mass. 
With the lowest order perturbation terms in fliOD , 

Q 2 v j 

Hi = 77 [(1 + In — )cos6 H sinO] , 

27reo r mcr 

where we have written down the vector products explicitly, we will now concentrate on 
the corrections in the energy, taking the harmonic oscillator excited states, (above the 
ground state), as the unperturbed ones. Using the notation of 0, let us consider only 
the ground state and the first excited state wave functions. 

a) Ground state: non-degenerate, E = HQ, n = 0, n r = 0, M = 0, 

fx Xr 2 

u Ofi (r,0) = \ -exp( — — ) , (48) 

V 7T L 

b) First excited state: doubly degenerate, E = 2hQ, n = 1, n r = 0, M = ±1, 

A Xr 2 

uo±(r,9) = —j=r exp( )exp(±i0) . (49) 

\/7T 2 



Here A = ^ has dimension of (length) 2 and n, n r and M are respectively the total, 
radial and orbital angular momentum quantum numbers, with M being an integer. 

It is readily seen that Hi above is linear in cosO and sinO and hence there are no first 
order effects in any of the states. 

Coming now to second order perturbation theory, we have 

. _ ^ |< m! I V I m >\ 2 

A 2 E m = Z m , ] - 1 _ ' { - . (50) 

For the ground state this reduces to 

_ l< 1 1 v 1 >l 2 l< -1 1 U I >| 2 

A2E ° ~ — w — w 1 w — w 

-C/0 hi l hlQ hl_ 1 
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Q2f] '" 2 [A 2 + -4^B 2 ] . (51) 



Eq — E\ Aneo m 2 c 2 
Here we have used the following relations, 

< m! | cosO I ' sinO | >= 5 m ',i ± ^m',-i , 

< ±1 I cosd I >= — / e- i{±e) cos9d9 = - , 

< ±1 I sinO I >= ±- . 
1 1 2 

A, (a dimensionless quantity), and B, (a quantity of dimension (length)' 1 ) , are defined 

as, 

A =< 1 I (1 + Zn— ) I >= 2(A) 3 /V(l + ln-)e- Xr2 dr , 
r Jo r 

B =< 1 I (-) I >= 2(A) i ( R re~ Xr2 dr . (52) 
r 7o 

We are coming to the value of R, the upper limit of integration. One can rewrite the 
result in ( |5lD as, 

- -ss0 2(Ar » m + ( ^ )2/22) • (53) 

where, 

Jj = f r ° k 2 (l + lnk)exp(-\r 2 k 2 )dk, (54) 

I 2 = f r ° ke(-\rlk 2 )dk = -^(1 - ea;p~^ . (55) 
Jo 2Arg 

The indefinite integral with the logarithm in ([54] ) is not obtainable in a closed form. We 
have the following identities 0, 

[* k 2 exp(-\r 2 k 2 )dk = i(Ar 2 )-f 7 (^, XR 2 ), (56) 

roc 1 , Q 

/ k 2 exp(-Xr 2 k 2 )dk = -(Ar 2 )^r(-), (57) 

/•OO A /7F S 

/ k 2 exp(-Xr 2 k 2 )ln(k)dk = ^— (Ar 2 )-^ (2 - Jn(4/i) - C), (58) 
jo 8 

where 7 and T are respectively the incomplete Gamma function and Gamma function 

and C is the Euler number. Note that in the limit R — ► 00, A2-E0 i n (|53|) becomes totally 

independent of ro. However, for any finite value of R, survives in the 7 function and 

in the integral in ([58]) as well. 

The non-trivial issue is the r dependence of the final result, which is coming from 

the upper limit ^. If R is infinity, then the r -dependence disappears but conceptually 

the behaviour of the particles is difficult to visualise: The attraction (repulsion) between 

opposite (similar) charges changes to repulsion (attraction), once the ro-barrier is crossed. 

It might seem natural to consider the upper bound of r to be r , but this compactification 

is also not satisfactory. 
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Generally one encounters logarithmic potentials in physical (3 dimensional) space hav- 
ing total infinite amount of charge with a translation symmetry in one space direction 
and the questions addressed here are not always pertinent. But if one is truely in a planar 
world, with logarithmic interparticle potential, then these questions can not be avoided. 
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Figure 3 
The functions, 

/(r) = y - ln(r) - J^, g(r) = Aee^r - |(1 + ^)][e* - r] 
are plotted for y = — 1, —2, —3, —4. 

(3a) y = -1 



(3b) y = -2 
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(3c) y 



(3d) y 
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Figure 1 

Potential profile in 3+1-dimensions, A = 0.01. 



Figure 2 

Potential profile in 2+1-dimensions, B = 10 -5 . 
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